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We show analytically and numerically that omission of quantum interference from the Chalker- 
Coddington model of the integer quantum Hall effect gives a localization length exponent v — 4/3 
as in ordinary two-dimensional percolation. Thus, contrary to semi-classical scaling arguments, 
tunneling alone does not lead to v = 7/3. 



Scaling and universality in the integer quantum Hall effect (IQHE) related to the localization-delocalization tran- 
sition of charge carriers when going from one plateau to the next is of considerable interest (l],|^] . By measuring how 
. ■ the width of the peaks of the longitudinal conductance scales with magnetic field B and linear size W of the system, 
Koch ct al. H found the scaling law W ~ \B — B c \~ v where v » 2.3 and B c is the value of the magnetic field needed 
to have half filling. 

This value excludes early theories |^,^| that predicted universal two-dimensional percolation exponent (y = i/ p = 
4/3). However, much before the measurements, Mil'nikov and Sokolov |J argued on the basis of a semi-classical 
picture and scaling that the classical value should be shifted by unity, i.e. v = 7/3. Numerically 2.35 was found 
in the Chalker-Coddington (CC) model 0, which includes in addition to the semi-classical tunneling effect also 
interference, and seems to capture the critical behavior in IQHE. It has been unclear whether tunneling is sufficient to 
explain, at least in a range of parameters, the change of the localization length exponent from the classical percolation 
value, 4/3, to the observed one, (~ 7/3) 

We present here a study of the effect of tunneling without interference in a CC type model. Our conclusions, based 
on analytic arguments and tested numerically, are as follows: 1) Tunneling alone does not change the localization 
-H , length exponent from its classical percolation value, 4/3. 2) The value actually observed is caused by the mechanism 
we left out (but emphasized by CC): Quantum interference. 

In the classical picture the derealization of the charge carriers is considered as a percolation transition 
Noninteracting charges are assumed to move in a disordered potential slowly varying on the scale of the magnetic 
length / = \/Ti/eB (high field model). In this approximation, the electron wave function tp is localized to equipotential 
curves, forming "ribbons" of width I. Changing B causes the localized wave function to move up or down in this 
landscape. Typically, the ribbons constituting the wave function are restricted to encircling mountains or valleys in 
the energy landscape. However, scattering between different ribbons is possible when they overlap at saddle points. 
There is a typical distance a between saddle points. In a window AB centered at B c , a sample-spanning cluster of 
ribbons overlapping at saddle points exists, resulting in derealization. The width of the window AB is determined by 
the magnetic length I and decreases with increasing sample size W. This classical description becomes semi-classical 
once the question of overlap of wave functions across saddle points takes tunneling into account. Overlap is then 
G , characterized by a tunneling probability fl2f 



p= 1 . ml 2 \V XX Vyy\ 1/2 

1 + eVx ' X he\B-B c \ ' 1 ' 



Here V xx and V yy are the principal double spatial derivatives of the disordered potential at the saddle point, and m is 
the effective electron mass. The semi-classical argument leading to v = v v + 1 = 7/3 is essentially based on a change 
of effective basic length scale of the percolation problem from a to x a - However, there is an underlying assumption 
that the tunneling does not change in any essential way the structure of the infinite cluster defined by the topology 
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of the paths in the random energy landscape. As we will show later, this is the weak part of the argument, and the 
cause of its breakdown. 

In order to proceed we study the semi-classical argument using the extended CC model Jl3| . We define first a fully 
classical (percolation) version |l4| ] . Imagine a square lattice of linear size L where every plaquette is occupied by one 
of the tiles of Fig. El The result is a set of lines which are either closed or end at the boundary of the sample, see 
Fig. ||a. A control parameter p tells the probability of the tile la. This problem can be mapped exactly to bond 
percolation on a square lattice fl4j . The mapping is accomplished by drawing a diagonal on each tile oriented in such 
a way that it does not intersect the two quarter circles already on the tile. The resulting pattern of diagonals is shown 
in Fig. ||a. This pattern consists of two lattices oriented 45° with respect to the original one defined by the tiles, and 
with lattice constant equal to the length of the diagonals of the tiles. The two lattices are displaced by a distance 
equal to the linear size of the tiles in the two principal orthogonal directions. Furthermore, whenever there is a bond 
present in one of the two lattices, the corresponding one is missing in the other lattice, and vice versa, making the 
two lattices duals of each other. This is shown in Figs. ||a, b and c. This makes the problem of random orienting 
the two tiles of Fig. |l| equivalent to placing a bond or not in one of the two intertwined lattices, say that of fig. ||b. 
It is also clear that, for p = 1/2 we deal with a percolation problem at the critical point, as with equal probability 
in choosing either tile corresponds to placing a bond with probability 1/2, which is the percolation threshold on the 
square lattice. 

The bonds of the tiling represent the classical paths the electrons may take in the disordered energy landscape. 
Each bond can only be traversed in one direction. The possible directions are chosen so that the network forms a grid 
of loops of alternating handedness according to the magnetic field and the surface topology. The centers of the the 
tiles represent saddle points, and the lattice constant corresponds to a. The question of delocalization is equivalent 
to whether there is a path spanning the system. 

We now introduce tunneling within the framework of a semi-classical approximation. The transmission probabilities 
at a given saddle point are controlled by a parameter 7 characterizing the amount of tunneling p3| . By adjusting it 
to the usual control parameter of the CC model determining the scattering properties at the saddle point we have 

(2) 
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When 7— ► or 7 — ► 00, tunneling is negligible, while when 7 = 7 C = ln(l + V2), the transmission probability is 1/2. 

Following Lee et al. fl3|] , we assume 7 = r y c e^~ v . Here p is the dimensionless chemical potential equal for all saddle 
points in the network, while v is chosen randomly on the interval [— w/2, +w/2] representing the random potential. 
In the case when \p — v\ <C 1, we may relate \p — v\ cx \B — B c \ of the high- field model. In extreme case when 
w — > 00, there is no contact across the saddle points and the model is just the classical version discussed above with 
1/2 occupation probability for both tiles. Delocalization happens then as a classical percolation phenomenon and, 
accordingly, the correlation length exponent has to be 4/3. 

We now turn to the case of finite w when tunneling becomes important. In the CC model one should assign phases 
to the electrons moving in the channels and the change of the phases during tunneling should be taken into account. 
We neglect the phases here since the aim of this work is to test the effect of tunneling alone. What is the action of 
tunneling near a given saddle point, i.e. on a given tile in this Fig. |^a? It simply means that, with some probability, a 
tiling in the classical picture has to be changed to the other one for this particular electronic process. This corresponds 
in the bond percolation picture either to remove or to place a bond with some well defined probability that mimics 
Eqs. (^). Thus, instead of introducing a 7 = 7 c exp(/i — v) and relating it to the tunnel probability P by Eq. (|J), we 
can divide the assigning construction of the conducting paths into two parts: one which places the "classical" paths, 
and then a second probability that models the tunneling. Given a control parameter p S [0, 1], we draw two random 
numbers r and p from a uniform distribution between zero and one for each bond. The bond is present if the following 
function is equal to one, 

7r(p, A; r, p) = 6(p - r) + (0(r - p) - 6(p - r)) 6{P - p) , (3) 

where, 9 is the step function and P is defined in Eq. ([!]), with 1/x = h\p — r| and A being a scale factor. In words: 
Tunneling introduces a bond if there was no one classically and cancels the bond if it was present. This means that 
actual realizations of paths are followed and an average over such paths should finally carried out. 

The smaller A, the more important is tunneling. The idea behind Eq. (|^) is simple. If A — > 00, P — > 0, and Eq. (|^) 
reduces to 7r(p; r, p) = 6(p — r), which is describing a simple bond percolation problem with p as control parameter. 
We may therefore interpret r as characterizing the level of the corresponding saddle point, and p the magnetic field. 
The second term in Eq. (|^) will switch the value of n given by the first "classical" term to the opposite value with a 
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probability which falls off exponentially the further p is away from the value r. In particular, when p = r for bond i, 
there is a 50% chance that the value of ir will be reversed compared to its classical value. Furthermore, the further r 
and p are apart, the more difficult to switch the value of ir. This completes the definition of our tunnel-only version 
of the Chalker-Coddington model. 

The probability that a bond is present is then 

1 , f 1 , , , 1 . ( 1 + e- A P 



?r(p) = (tt(>; r, p)) = dr dp vr(p; r, p) = p + - log I - - e _ A(1 _ rt J , (4) 

Since assigning a bond or not is a purely local decision, which is so far based on two random numbers r and p, we 
may simplify and draw a single random number q between zero and one and comparing this to W: 

wfaq)=0(W-q). (5) 

This formulation is equivalent to that of Eq. (|^). ft is simpler to analyze, but harder to interpret. The effect of 
tunneling has been to produce a non-linear relationship between the control parameter p and the probability for a 
bond to be present. However, we have that W = 1/2 when p = 1/2, as expected: The system has the same critical 
point with or without tunneling: p c = 1/2. 

We are now ready to analyze what happens close to the critical value of p, and in particular determining the 
localization length exponent v in this model. In order to do this, we first note that using 7f as control parameter and 
not p, makes the problem a standard percolation problem. Thus, the correlation length £ diverges as 

^ ~ ]W C — tt]-"- , (6) 
where W c = 1/2, and v v = 4/3. We now expand W in powers of Sp — p c — p to get 

n(l/2 + Sp) = 1/2 + tanh Sp + 0{5p 3 ) . (7) 

Plugging this expansion into Eq. (0) gives 

t~\Pc-p\-">, (8) 
but with a different prefactor. The exponent in this equation is by definition v and we have shown that 
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v = v v = - (9) 

in the CC model when only tunneling is taken into account. Based on this result, we conclude that the semi-classical 
argument giving v = v v + 1 fails. 

It is interesting to note that in a one-dimensional version of this model, first studied in |]ll"|| , there is a shift of v in 
comparison to v p = 1 by precisely one fl5f| . Thus, in this case, the semiclassical argument does work. The difference 
between the one and two-dimensional cases is that in the latter the number of bonds almost belonging to the infinite 
cluster is so large — their density is finite — that the system simply is driven away from criticality. 

We have carried out numerical simulations of the semi-classical CC model as defined above in order to see the way 
asymptotics sets in. We generated different systems sizes up to 1100 x 1100 with at least 20 000 samples using Monte 
Carlo renormalization group fl6|| to determine the exponents. We have based our simulations on Eq. (|^) directly, 
and on Eq. (^) directly. In the latter case, we used geometrical sample averaging in order to test whether the typical 
percolation probability behaves differently from the average one. We found no difference in the critical behavior. The 
measured effective exponent was v — 1.53 for W = 100 and v = 1.45 for W = 1100. Thus there is a slow convergence 
towards the predicted universal exponent 4/3 but even small size samples have an effective exponent much smaller 
than 7/3, confirming our analytical conclusion. 

The localization length exponent observed in the Chalker-Coddington model and in other models of the integer 
quantum Hall effect in addition to experiments, is larger than v p by about unity. What we have shown here is that 
the tunneling mechanism alone is not able to account for the change of the exponent in any range of the parameters 
- therefore interference must be responsible for it. 

As the control parameter AB increases the interference becomes less and less pronounced and the entangled paths 
due to tunneling percolation dominate. This should be reflected in a crossover from the quantum exponent v w 2.3 
to the classical value 4/3. We suggest the experimental check of this predition. 
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FIG. 1. The two fundamental tiles constituting the building blocks of the CC model in the classical limit. 

FIG. 2. (a) A set of dragon curves with corresponding diagonals inscribed for p — 1/2. In (b) and (c) we have decomposed 
the set of diagonals into a sublattice (b) and the corresponding dual lattice (c). 
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